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Abstract 

A model-independent search for deviations from the Standai^d Model prediction is per- 
formed in e'^p and e~p collisions at HERA using HI data corresponding to an integrated 
luminosity of 117 pb~^. For the first time all event topologies involving isolated elec- 
trons, photons, muons, neutrinos and jets with high transverse momenta are investigated 
in a single analysis. Events ai^e assigned to exclusive classes according to their final state. 
A statistical algorithm is developed to search for deviations from the Standard Model in 
the distributions of the scalar sum of transverse momenta or invariant mass of final state 
particles and to quantify their significance. A good agreement with the Standard Model 
prediction is observed in most of the event classes. The most significant deviation is found 
for a topology containing an isolated muon, missing transverse momentum and a jet, con- 
sistent with a previously reported observation. 
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1 Introduction 



At HERA electrons' and protons collide at a centre-of-mass energy of up to 319 GeV. These 
high-energy electron-proton interactions provide a testing ground for the Standard Model (SM) 
complementary to e+e" and pp scattering. It is widely believed that the SM is incomplete and 
that new physics signals may appear below energies of 1 TeV. Many extensions to the SM have 
been constructed during the last decades predicting various phenomena which may be visible 
at high energies or large transverse momenta (Pt)- HERA data have been used to test some of 
these models of new processes by analysing their anticipated experimental signatures and limits 
on their parameters have been derived [1]. 

The approach described in this paper consists of a comprehensive and generic search for de- 
viations from the SM prediction at large transverse momenta. All high Pt final state configura- 
tions involving electrons (e), muons (fi), jets (j), photons (7) or neutrinos (v) are systematically 
investigated. The analysis covers phase space regions where the SM prediction is sufficiently 
precise to detect anomalies and does not rely on assumptions concerning the characteristics of 
any SM extension. Such a model-independent approach might discover unexpected manifesta- 
tions of new physics. Therefore it addresses the important question of whether evidence for new 
physics might still be hidden in the data recorded at collider experiments. A similar strategy for 
a model-independent search was previously presented in [2]. 

All final states containing at least two objects (e, /i, j, 7, i^) with Pt > 20 GeV in the polar 
angle^ range 10° < 6* < 140° are investigated. The complete HERA I data sample (1994 - 2000) 
is used, corresponding to an integrated luminosity of 1 17 pb"'^. All selected events are classified 
into exclusive event classes according to the number and types of objects detected in the final 
state (e.g. e-j, jJi-j-y, j-j-j-j-j)- These exclusive event classes ensure a clear separation of final 
states and allow an unambiguous statistical interpretation of deviations. All experimentally 
accessible combinations of objects have been studied and data events are found in 22 of them. 

In a first analysis step the global event yields of the event classes are compared with the SM 
expectation. The distributions of the invariant mass Mau and of the scalar sum of transverse 
momenta ^ Pt of high Pt final state objects are presented. New physics may be visible as an 
excess or a deficit in one of these distributions. Therefore, in a second step these distributions 
are systematically investigated using a dedicated algorithm which locates the region with the 
largest deviation of the data from the SM prediction. The probability of occurrence of such a 
deviation is derived, both for each event class individually and globally for all classes combined. 

This paper is organised as follows. Section 2 describes the Standard Model processes at 
HERA and their Monte Carlo simulation. The HI detector, the event selection and measurement 
procedure are described in section 3. The event yields and distributions for each event class are 
presented in section 4. The search strategy and results are explained in section 5. Section 6 
summarises the paper. 

'in this paper "electrons" refers to both electrons and positrons, if not otherwise stated. 

^The origin of the HI coordinate system is the nominal ep interaction point, with the direction of the proton 
beam defining the positive 2-axis (forward region). The transverse momenta are measured in the xy plane. The 
pseudorapidity rj is related to the polar angle dhy i] = — In tan(0/2). 
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2 Standard Model processes and Monte Carlo generation 



Several Monte Carlo event generators are combined to simulate events for all SM processes 
which have large cross sections or are expected to be dominant for specific event classes, avoid- 
ing double-counting. All processes are generated with an integrated luminosity significantly 
higher than that of the data sample and events are passed through a full detector simulation [3]. 
At high transverse momenta the dominant SM processes are the photoproduction of two jets 
and neutral current (NC) deep-inelastic scattering (DIS). In the following the abbreviation X 
represents all reaction products other than the high Pt objects considered. 

Photoproduction of jets and photons To simulate the direct and resolved photoproduction 
of jets ep — > jjX, prompt photon production ep — » 'yjX and the resolved photoproduction 
of photon pairs ep 77X, the PYTHIA 6.1 event generator [4] is used. Light and heavy 
flavoured jets are generated. The simulation contains the Bom level hard scattering matrix 
elements and radiative QED corrections. 

Neutral current deep-inelastic scattering The Bom, QCD Compton and Boson Gluon Fu- 
sion matrix elements are used in the RAPGAP [5] event generator to model NC DIS events. The 
QED radiative effects arising from real photon emission from both the incoming and the outgo- 
ing electrons are simulated using the HERACLES [6] generator. Hence the NC DIS prediction 
contains the processes ep ejX, ep — > ejjX and also models final states with an additional 
radiated photon. 

Charged current deep-inelastic scattering Charged current (CC) DIS events are simulated 
using the DJANGO [7] program, which includes first order QED radiative corrections based on 
HERACLES. This prediction contributes to the final states ep — > I'jX, ep — ^ i^jjX and to final 
states with an additional radiated photon. 

QED Compton scattering Elastic and quasi-elastic Compton processes ep e'jX are sim- 
ulated with the WABGEN [9] generator. The inelastic contribution is already included in the 
NC DIS RAPGAP sample. 

Electroweak production of lepton pairs Multi-lepton events (ee, fifi, tt) are generated with 
the GRAPE [10] program, which includes all electroweak matrix elements at tree level. Multi- 
lepton production via 77, 7Z, ZZ collisions, intemal photon conversion and the decay of virtual 
or real Z bosons is considered. Initial and final state QED radiation is included. The complete 
hadronic final state is obtained via interfaces to PYTHIA and SOPHIA [11] for the inelastic and 
quasi-elastic regimes, respectively. 

W production The production of W bosons ep — > WX and ep WjX is modelled using 
EPVEC [12]. Next-to-leading order QCD corrections [13] are taken into account by reweighting 
the events as a function of the transverse momentum and rapidity of the W boson [14]. 
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Processes with the production of three or more jets, e.g. ep — > iijX or ep —* jjjjX, are ac- 
counted for using leading logarithmic parton showers as a representation of higher order QCD 
radiation, with the exception of CC DIS, where the colour-dipole model [8] is used. Hadronisa- 
tion is modelled using Lund string fragmentation [4]. The prediction of processes with two or 
more high transverse momentum jets, e.g. ep jjX, ep — > ejjX, is scaled by a factor of 1.2 
to normalise the leading order Monte Carlos to next-to-leading order QCD calculations [15]. 



3 Experimental technique 

3.1 The HI detector 

The HI detector [16] components relevant to the present analysis are briefly described here. 
Jets, photons and electrons are measured with the Liquid Argon (LAr) calorimeter [17], which 
covers the polar angle range 4° < 9 < 154° with full azimuthal acceptance. Electromag- 
netic shower energies are measured with a precision of aiE)/E = 12%/^E/GeY © 1% and 
hadronic energies with a{E)/E = 50%/^E/GeV © 2%, as measured in test beams. The 
central and forward tracking detectors are used to measure charged particle trajectories, to re- 
construct the interaction vertex and to supplement the measurement of the hadronic energy. 
The innermost proportional chamber CIP (9° < 9 < 171°) is used to veto charged particles 
for the identification of photons. The LAr and inner tracking detectors are enclosed in a super- 
conducting magnetic coil with a strength of 1.15 T. The return yoke of the coil is the outermost 
part of the detector and is equipped with streamer tubes forming the central muon detector 
(4° < 6* < 171°). It is also used to supplement the measurement of hadrons. In the forward 
region of the detector (3° < 6^ < 17°) a set of drift chamber layers (the forward muon system) 
detects muons and, together with an iron toroidal magnet, allows a momentum measurement. 
The luminosity measurement is based on the Bethe-Heitler process ep ep'j, where the photon 
is detected in a calorimeter located downstream of the interaction point. 

The main trigger for events with high transverse momentum is provided by the LAr calorime- 
ter. The trigger efficiency is close to 100% for events having an electromagnetic deposit in the 
LAr (electron or photon) with transverse momentum greater than 20 GeV [19]. Events trig- 
gered only by jets have a trigger efficiency above 90% for P^^'^ > 20 GeV and nearly 100% 
for P^^^ > 25 GeV [20]. For events with missing transverse momentum above 20 GeV, deter- 
mined from an imbalance in the transverse momentum measured in the calorimeter, the trigger 
efficiency is ~ 90%. The muon trigger is based on single muon signatures from the central 
muon detector, which are combined with signals from the central tracking detector. The trigger 
efficiency for di-muon events is about 70% [21]. 

3.2 Event selection 

At HERA electrons or positrons with an energy of 27.6 GeV collide with protons at an energy of 
920 GeV resulting in a centre-of-mass energy of y/s = 319 GeV. Before 1998 the proton energy 
was 820 GeV resulting in a centre-of-mass energy of ^/s = 301 GeV. The event sample studied 
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consists of the full 1994-2000 HERA I data set. It corresponds to an integrated luminosity of 
36.4 pb~^ in e^p scattering at -^/i = 301 GeV and 13.8 pb"^ in e~p scattering and 66.4 pb~^ in 
e+p scattering at = 319 GeV. 

The data selection requires at least one isolated electromagnetic cluster, jet or muon to be 
found in the detector acceptance. Energy deposits in the calorimeters and tracks in the in- 
ner tracking system are used to form combined cluster-track objects, from which the hadronic 
energy is reconstructed. To reduce background it is demanded that the event vertex be re- 
constructed within 35 cm in z of the nominal position^ and that {Ei — P^^i) < 75 GeV, 
where Ei is the particle energy and Pz,i is the z component of the particle momentum. Here, 
the index i runs over all hadronic energy deposits, electromagnetic clusters and muons found 
in the event. Due to energy-momentum conservation events are expected to have a value of 

{Ei — Pz^i) = 55.2 GeV, twice the electron beam energy, if only longitudinal momentum 
along the proton beam direction is unmeasured. Events with topologies typical of cosmic ray 
and beam-induced background are rejected [22]. Moreover, the timing of the event is required 
to coincide with that of the ep bunch crossing. 

The identification criteria for each type of particle are based on those applied in previous 
analyses of specific final states [15, 19,21,23]. Additional requirements are chosen to ensure 
an unambiguous identification of particles, whilst retaining high efficiencies. The following 
paragraphs describe the identification criteria for the different objects and give the identification 
efficiencies for the kinematic region considered in the analysis. 

Electron identification The electron identification is based on the measurement of a com- 
pact and isolated electromagnetic shower in the LAr calorimeter. The hadronic energy within 
a distance in the pseudorapidity-azimuth (rj — 0) plane R = (A?])^ + (A0)2 < 0.75 around 
the electron is required to be below 2.5% of the electron energy. This calorimetric electron 
identification is complemented by tracking conditions. A high quality track is required to geo- 
metrically match the electromagnetic cluster within a distance of closest approach to the cluster 
centre-of-gravity of 12 cm. No other good track is allowed within i? < 0.5 around the electron 
direction. In the central region (20° < 9 < 140°) the distance between the first measured point 
in the central drift chambers and the beam axis is required to be below 30 cm in order to reject 
photons that convert late in the central tracker material. In addition, the transverse momentum 
measured from the associated track P^^'' is required to match the calorimetric measurement P^ 
with 1/P^*'' — l/Px < 0.02 GeV~^. In the region not fully covered by the central drift cham- 
bers (10° < 9 < 37°) a wider isolation cone of _R = 1 is required to reduce the contribution 
of fake electrons from hadrons. In this forward region the identification is completed by the 
requirement of associated hits in the CIP, which reduces the contamination from neutral parti- 
cles showering in the material of the forward region. The resulting electron finding efficiency 
is 85% in the central region and 70% in the forward region. 

Photon identification The photon identification relies on the measurement of an electromag- 
netic shower and on the same calorimetric isolation criteria against hadrons as for the electron 
identification. In addition, photons are required to be separated from jets with P;jf^ > 5 GeV 

^This is not required for the event classes containing only photons or photons and a neutrino. 
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by a distance of i? > 1 to the jet axis. Vetoes are applied on any charged track pointing to the 
electromagnetic cluster. No track should be present with a distance of closest approach to the 
cluster below 24 cm or within R < 0.5. An additional veto on any hits in the CIP is applied. 
The resulting photon identification efficiency as derived using elastic QED Compton events is 
90%. 

Muon identification The muon identification is based on a track in the forward muon system 
or in the inner tracking systems associated with a track segment or an energy deposit in the 
central muon detector [23]. The muon momentum is measured from the track curvature in the 
toroidal or solenoidal magnetic fields. A muon candidate should have no more than 8 GeV 
deposited in the LAr calorimeter in a cylinder of radius 0.5 in (77 — 0) space, centred on the 
muon track direction. In di-muon events, the requirement of an opening angle between the 
two muons smaller than 165° discards cosmic ray background. Beam halo muons are rejected 
by requiring that the muons originate from the event vertex. Finally, misidentified hadrons are 
almost completely suppressed by requiring that the muon candidate is separated from the closest 
jet with P^,"^* > 5 GeV by i? > 1. The efficiency to identify muons is greater than 90% [23]. 

Jet identification Jets are defined using the inclusive fc^ algorithm [24,25]. The algorithm 
is applied in the laboratory frame with a separation parameter of 1 and using a Pt weighted 
recombination scheme [24] in which the jets are treated as massless. The jet algorithm is run 
on all combined cluster-track objects not previously identified as electron or photon candidates. 
The scattered electron may fake a jet. This effect is important for multi-jet events, especially at 
high transverse momenta. To reject these fake jets, the first radial moment of the jet transverse 
energy [26,27] is required to be greater than 0.02 and the quantity M^^^ /P^^ must be greater 
than 0.1 [15,27]. The invariant mass M^^^ is obtained using the four-vector of all objects 
belonging to the jet. If the fraction of the jet energy contained in the electromagnetic part 
of the LAr calorimeter is greater than 0.9, the above criteria are tightened to 0.04 and 0.15, 
respectively. The jet selection efficiency is 97%. 

Neutrino identification A neutrino candidate is defined in events with missing transverse 
momentum above 20 GeV. The missing momentum is derived from all identified particles and 
energy deposits in the event. Fake missing transverse momentum may also arise from the 
mismeasurement of an identified object. This effect is reduced by requiring that the neutrino^ 
be isolated from all identified objects with a transverse momentum above 20 GeV. Requiring 
{Ei — Pz^i) < 48 GeV discards neutrino candidates from NC processes where the missing 
transverse momentum is caused by energy leakage in the forward region. If exactly one electron 
or muon object is found, a neutrino object is only assigned to an event if A0(/ — Xtot) < 170°, 
where A0(/ — Xtot) is the separation in azimuthal angle between the lepton / and the direction 
of the system Xtot built of all hadronic energies. 

^The four- vector of the neutrino is calculated under the assumption {Ei — Pz,i) + {Ei, — Pz,v) — 55.2GeV. 
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Event classification The common phase space for electrons, photons, muons and jets is de- 
fined by 10° < 9 < 140° and > 20 GeV. The neutrino phase space is defined by missing 
transverse momentum above 20 GeV and ^ . [Ei — Pz,i) < 48 GeV. These values are chosen 
to retain a high selection and trigger efficiency. All particles with Pt > 20 GeV, including 
the neutrino defined by its reconstructed four-vector, are required to be isolated compared with 
each other by a minimum distance R of one unit in the r/ — plane. The events are classified, 
depending on the number and types of objects, into exclusive event classes. Events with an 
isolated calorimetric object in the considered phase space which is not identified as a photon, 
electron or jet are discarded from the analysis in order to minimise wrong classifications. 

Based on these identification criteria, purities can be derived for each event class with a 
sizeable SM expectation. Purity is defined as the ratio of SM events reconstructed in the event 
class in which they are generated to the total number of reconstructed events in this class. Most 
purities are found to be above 60% and they are close to 100% for the j-j, e-j, j-u and event 
classes. 



3.3 Systematic uncertainties 

This section describes the sources of experimental and theoretical systematic uncertainties con- 
sidered. Experimental systematic uncertainties arising from the measurement of the objects are 
listed in table[l](for more details see [20, 27, 28]). 

• The electromagnetic energy scale uncertainty varies between 0.7% and 3% depending 
on the particle's impact point on the LAr calorimeter surface [19]. The polar angular 
measurement uncertainty of electromagnetic clusters varies depending on 9 between 1 
and 3 mrad [19]. The identification of electron and photon candidates depends on the 
tracking efficiency, which is known with a precision ranging from 2% for polar angles 
above 37° to 7% in the forward region. 

• The hadronic energy scale of the LAr calorimeter is known to 2%. The uncertainty on the 
jet polar angle determination is 5 mrad for 9 < 30° and 10 mrad for 9 > 30°. 

• The uncertainty on the transverse momentum measurement for muons amounts to 5%. 
The uncertainty on the polar angle is 3 mrad. The muon identification efficiency is known 
with a precision of 5%. 

• The trigger uncertainties for each class are determined by the object with the highest 
trigger efficiency. The uncertainty on the trigger efficiency is estimated to be 3% if the 
event is triggered by a jet or missing transverse momentum and 5% if it is triggered by a 
muon. For electrons and photons the uncertainty on the trigger efficiency is negligible. 

• The uncertainty in the integrated luminosity results in an overall normalisation error of 
1.5%. 

• The uncertainty in the reconstruction of J2i i^i " Pz,i) and the missing P^ for the neu- 
trino classification are obtained by propagation of the systematic errors for other objects. 
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Depending on the generator level production process, different theoretical uncertainties are 
used as listed in table |21 The errors attributed to the predictions for ep jjX, ep j'jX, 
ep — ^ ji^X, ep jeX, ep — > jji^X, ep jjeX and W production include uncertainties in the 
parton distribution functions and those due to missing higher order corrections [15, 23, 27, 28]. 
The error attributed to ep jijjiX and ep eeX results mainly from the limited knowledge 
of the proton structure [21, 29]. The error on the QED Compton cross section is estimated to 
be 5% for elastic and 10% for inelastic production. An additional theoretical error of 20% is 
applied for each jet produced by parton shower processes (e.g. 20% for the j-j-j event class). 
An uncertainty of 50% is added to the prediction for NC DIS events with missing transverse 
momentum above 20 GeV and a high electron. This uncertainty is estimated by a compar- 
ison of the missing transverse momentum distribution between NC DIS events with a low Pt 
electron {Pt < 20 GeV) and the SM prediction. 

All systematic errors are added in quadrature and are assigned to the SM predictions. For 
example, the resulting total uncertainties for e-j events are 10% and 35% at low and high in- 
variant mass Mall, respectively. In the j-j event class the errors are typically 20% and reach 
40% - 50% for Mall and ^ Pt values around 250 GeV. 



4 Event yields 

All possible event classes with at least two objects are investigated^. The event yields subdivided 
into event classes are presented for the data and SM expectation in figure [T] All event classes 
with a SM expectation greater than 0.01 events are shown. No other event class contains data 
events. The distributions of the scalar sum of transverse momenta ^ Pt and of the invariant 
mass Mall of all objects are presented in figures|2land|3]for classes with at least one event. 

The dominant high Pt processes at HERA, i.e. photoproduction of jets, NC and CC DIS, 
dominate in the j-j, e-j and j-u event classes, respectively. Events are observed with ^Pt 
and Mall values as large as 250 GeV. A good description of the data spectra by the prediction is 
observed. The prediction for the event classes j-j-7 and e-j-7 is dominated by photoproduction 
and NC DIS processes with the radiation of a photon, respectively. There is good agreement 
between the data and the prediction. No event is observed in the radiative CC classes z/-7 and 
j-z/-7, consistent with the expectation of 2.1 ± 0.3 and 1.0 ± 0.1, respectively. The j-j-j, e-j-j, 
e-j- j-j, j-j-y and j-j-j-y event classes correspond to processes with additional jet production 
due to higher order QCD radiation. The yields of these event classes are also well described by 
the SM prediction. 

The e-7 event class is dominated by QED Compton scattering processes and ^ Pt and Mau 
values up to 160 GeV are measured. A good agreement with the SM is observed. The prompt 
photon j-7 event class extends up to Mau ~ 150 GeV and is well described by the prediction. 
The purity in this class is moderate (40 — 50%) due to the high background from misidentified 
electrons in NC DIS. Backgrounds where hadrons are misidentified as photons are small. One 

^The fi-v event class is discarded from the present analysis. It is dominated by low transverse energy photo- 
production events in which a poorly reconstructed muon gives rise to missing transverse momentum, which fakes 
the neutrino signature. 
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event is observed in the 7-7 event class for an expectation of 1.1 ± 0.5, which is dominated by 
the ep — > e'yX process, where the electron is misidentified. Contributions of higher order QED 
processes, which could lead to two high transverse momentum photons, are negligible. 

Lepton pair production from 77 processes dominates in event classes with several leptons. 
The e-e event class contains 8 events for an expectation of 11.2 ± 1.4. In this channel, a dis- 
crepancy with the SM expectation was previously reported for high masses by the HI collabo- 
ration [29]. All multi-electron events mentioned in [29] and located in the phase space of this 
analysis are found. In the region Mau > 100 GeV, 3 events are observed and 1.16 ± 0.25 are 
expected. The higher SM prediction compared with the prediction of 0.3 in [29] is due to the 
extended polar angle range in the forward region. This leads to an additional ^0.4 background 
events from fake electrons and ^ 0.4 events from the ep eeX processes. The e-e-e class con- 
tains no events. None of the tri-electron events of [29] are selected here due to the requirement 
of high transverse momentum. The predictions for the e-/i and event classes are dominated 
by muon pair production from two-photon reactions. The e-fi event class is populated if the 
scattered electron and only one of the muons are selected. In the e-/i class, 4 events are ob- 
served compared with an expectation of 4.8 ± 0.6. A slight excess is observed in the fi-fi event 
class where 6 events are found and 2.7 ± 0.6 are expected. Muon pair production processes also 
contribute ~ 85% in the /i-j event class, where a good agreement is found. In the e-/i, fi-fi and 
yu-j event classes the ^ Pt and Mau values of the data lie between 50 and 100 GeV. 

The prediction for the event classes fi-j-i' and e-j-u consists mainly of high Pt W produc- 
tion with a subsequent leptonic decay. A discrepancy between the data and the SM expectation 
is observed in the fi-j-i' event class, where 4 events are observed for an expectation of 0.8 ± 0.2. 
The Yl Pt values reach 170 GeV and the Mau values 200 GeV. In this event class less than 
0.002 background events are expected from the photoproduction of jets via QCD processes. 
Such a deviation was previously reported in [23] and will be further discussed in Section 5. In 
the e-j-u event class 2 data events are observed for an expectation of 0.9 ± 0.2. Some of the 
e-j-u events mentioned in [23] have an electron with a transverse momentum below 20 GeV and 
are therefore not selected as e-j-u events in the present analysis. The event topology e-u is also 
expected to contain events arising from W production together with background from NC DIS. 
In the e-u event class, 9 data events are observed compared with an expectation of 12.9 ± 4.5. 

A slight excess of the data compared with the prediction is observed in the j-j-j-j event 
class, with 10 data events observed and 5.2 ± 2.2 expected. One event is observed in the 
e- j-j-j-j event class, to be compared with an expectation of 0.026 ± 0.011. This event has 
a of 207 GeV and an invariant mass Maii of 262 GeV. The NC DIS expectation for 

Mall > 260 GeV is (9 ± 6) ■ 10"^ as derived using RAPGAP. The energy flow of the event in 
the r] — (p view is presented in figure HI The NC DIS and photoproduction SM predictions have 
been tested using a sample of j-j-j-j events with P^*^* > 15 GeV and e-j-j-j-j events with 

> 10 GeV and P|>^* > 5 GeV. An adequate description of the ^ Pt and Mau distributions 
of the data is obtained within the quoted SM uncertainties. Since the NC DIS prediction for 
Mall > 260 GeV is only of order 0.001 fb, rare SM processes not considered in this analysis 
such as W pair production may be dominant in this kinematic domain. 

No events are found in any other event class, in agreement with the SM expectation (see 
figure [T]). 
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5 Search for deviations from the Standard Model 
5.1 Search algorithm and strategy 

In order to quantify the level of agreement between the data and the SM expectation and to 
identify regions of possible deviations, a new search algorithm is developed. Detailed studies 
have shown that Man and E Pt have a large sensitivity to new physics (see appendix and [28]). 
The algorithm described in the following locates the region of largest deviation of the data from 
the SM in these distributions. The calculation of the significance of this deviation is inspired 



Definition of regions A region in the E -Pr and Mau distributions is defined as a set of con- 
nected histogram bins^ with a size of at least twice the resolution. All possible regions of any 
width and at any position in the histograms are considered. The number of data events (Nobs), 
the SM expectation (Nsm) and its total systematic uncertainty (5Nsm) are calculated for each 
region. 

Determination of the most interesting region A statistical estimator p is defined to judge 
which region is of most interest. This estimator is derived from the convolution of the Pois- 
son probability density function (pdf) to account for statistical errors with a Gaussian pdf, 
G{b; Nsm, ^Nsm), with mean Nsm and width SNsm, to include the effect of non negligible 
systematic uncertainties. The estimator is defined via 



The factor A ensures normalisation to unity. If G is replaced by a Dirac delta function 6{b — 
Nsm) the estimator p becomes the usual Poisson probability. The value of p gives an estimate of 
the probability of a fluctuation of the SM expectation upwards (downwards) to at least (at most) 
the observed number of data events in the region considered. The region of greatest deviation is 
the region having the smallest p-value, Pmin- Such a method is able to find narrow resonances 
and single outstanding events as well as signals spread over large regions of phase space in 
distributions of any shape [28]. 

*In order to minimise binning effects, a bin size smaller than or comparable with the resolution of the studied 
quantity is used. A 5 GeV bin size is used for all distributions. Further reduction of the bin size has a negligible 
effect on the results. 



by [2]. 
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Significance per event class The probability that a fluctuation with a p- value at least as small 
as Pijyji occurs anywhere in the distribution is estimated using the following method. Many 
independent hypothetical data histograms are made by filling each bin with an event number 
diced according to the pdfs of the SM expectation (again a convolution of Poisson and Gaussian 
pdfs). For each of these hypothetical data histograms the algorithm is run to find the region of 
greatest deviation and the corresponding is calculated. The probability P is then defined 
as the fraction of hypothetical data histograms with a p^^^ equal to or smaller than the p^nin 
value obtained from the data. P is a measure of the statistical significance of the deviation 
observed in the data. If the event classes are exclusive, the P values can be used to compare the 
results of different event classes. Depending on the final state, a pj^in-value of 5.7 ■ 10^ ("5ct") 
corresponds to a value of — logio ^ between 5 and 6. 

Global significance The overall degree of agreement with the SM can be further quantified 
by taking into account the large number of event classes studied in this analysis. The probability 
of observing an event class with a given P value or smaller can be calculated with Monte Carlo 
(MC) experiments. A MC experiment is defined as a set of hypothetical data histograms (either 
in Mall or in J2 Pt) following the SM expectation with an integrated luminosity of 1 17 pb~^, 
on which the complete search algorithm and statistical analysis are applied as for data. This 
procedure is repeated many times. The expectation for the P values observed in the data is then 
given by the distribution of P^^^ values obtained from all MC experiments. The probability 
to find a P value smaller than the minimum observed in the data can thus be calculated and 
quantifies the global significance of the observed deviation. 

5.2 Search results 

The final P values obtained for event classes having at least one observed event are summarised 
in tabled The regions selected by the algorithm are presented in figuresEland|3J 

The most significant deviation of the analysis is found in the n-j-iy event class. This class 
has P values of 0.010 (Mau) and 0.001 (^ Pt). The mass region (155 < Mau < 200 GeV) 
contains 3 data events for an expectation of 0.19 ± 0.05. In the chosen ^ Pt region (145 < 
J2 Pt < 170 GeV) three data events are found while only 0.07 ± 0.03 are expected. This event 
topology was studied in [23] where this deviation at high Pt was already reported. 

A P value of 0.019 is found in the e-e event class in a region at high transverse momenta, 
100 < J2 < 130 GeV where 3 events are observed for an expectation of 0.18 ± 0.08. The 
deviation is less prominent in the region selected in the invariant mass distribution due to a 
higher background from NC DIS events. This corresponds to the excess of data events also 
identified in [29] . 

A deficit is observed in the e-j event class in the Yl Pt distribution in the region 180 < 
Y^Pt < 210 GeV. For a SM expectation of 31.2 ± 5.0 only 12 data events are observed. The 
derived P value is 0.021. 

Due to the uncertainties in the SM prediction in the j-j-j-j and e-j-j-j-j event classes at 
the highest Maii and ^ Pt, where data events are observed (see section 4), no reliable P values 
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can be calculated for these classes. Consequently, these event classes are not taken into account 
to determine the overall degree of agreement between the data and the SM. 

The P values for event classes with no data event observed and a small SM expectation 
are 1. This remains the case if an additional contribution is added from a possible further rare 
process not included here. Such classes can thus be considered in the calculation of the global 
significance. 

The P values observed in the data in all event classes are compared in figure |5l with the 
distribution of P^^^ obtained from the large set of MC experiments, normalised to one experi- 
ment. The comparison is presented for the scans of the Man and ^ Pt distributions. Most P 
values lie above 0.01, corresponding to event classes where no significant discrepancy between 
the data and the SM expectation is observed. The global probabilities to find at least one class 
with a P value smaller than the observation in the fi-j-u channel are 3% and 28% for the ^ Pt 
and Mall distributions, respectively (see appendix for details). 

To test the dependence of the analysis on the a priori defined Pt cuts, the whole analysis 
is repeated with two other object Pt cuts. The Pt cut was raised to 40 GeV for all objects 
and lowered to 15 GeV. In the latter case it was still required that at least one object has a 
Pt larger than 20 GeV in order to maintain a high trigger efficiency. The analysis was also 
repeated separately on the e^p and e^p data samples. In these four test scenarios a similar 
overall agreement with the SM is observed. The /i-j-i^ event class remains the one with the 
smallest P value in the scenario with a lowered Pt cut in the e^p data sample and no new 
discrepancy is observed. When raising the Pt cut to 40 GeV, it is mainly the two particle 
event classes containing jets that are still populated and the largest deviation is observed in the 
e-e class with P = 0.01. 



6 Conclusions 



The data collected with the HI experiment during the years 1994-2000 (HERA I) have been in- 
vestigated in a search for deviations from the SM prediction at high transverse momentum. For 
the first time all event topologies involving isolated electrons, photons, muons, neutrinos and 
jets are investigated in a single analysis. A good agreement between the data and the SM ex- 
pectation is found in most event classes. A better knowledge of rare processes may be required 
to search for deviations from the SM in final states with four jets at the highest invariant mass 
or transverse momentum. The distributions in the invariant mass and scalar sum of transverse 
momenta of the particles in each event class have been systematically searched for deviations 
using a statistical algorithm. The most significant deviation is found in the /i-j-i^ event class, a 
topology where deviations have also been previously reported. About 3% (28%) of hypothetical 
Monte Carlo experiments would produce a deviation in at least one event class which is more 
significant than that observed in the corresponding sum of transverse momenta (invariant mass) 
distribution of the topology with a jet, a muon and a neutrino. 
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Appendix 

Signals for new physics may appear either as a single deviation or a small set of deviations. 
The following outlines how a significant deviation might be defined and presents tests of the 
sensitivity of this analysis to specific signals for new physics. 

The probability to observe in the data a — log^^Q P greater than X in at least n event 
classes is given by the fraction of MC experiments having at least n event classes with a 

— logio -P"^^^ > X. The Px values obtained for this analysis are presented in tabled Up to 3 
event classes are considered. Since very similar P^ values are found for the Man and ^ P^ dis- 
tributions, averaged values are presented. For example, a P^ value smaller than 0.0005, which 
might be considered to represent a significant deviation, could be obtained from one event class 
with a — logj^Q P > 5, two event classes with a — log^^Q P > 3.5 or three event classes with a 

— logio -P > 3. It was found that one of these cases occurs either in Mau or ^ Py in around 
0.1% of all MC experiments. 

A set of pseudo data samples has been produced to test the sensitivity of the analysis pro- 
cedure to some dedicated signals for new physics. The prediction of a specific model for new 
physics is added to the SM prediction and this new total prediction is used to generate pseudo 
data samples. Again a Monte Carlo technique is used to vary the distribution of signal events 
and generate many MC experiments. The complete algorithm is run on those MC experiments 
and the mean value of — log^o P in all of them is derived as a measure of sensitivity of this 
analysis. 

The exotic production of top quarks via a flavour-changing neutral current is first investi- 
gated. The decay t bW with subsequent leptonic and hadronic W decays has been consid- 
ered. The (— log^Q P) values obtained are displayed in figure |3 (top) as a function of the cross 
section for producing a top when the proton beam energy is 920 GeV. Whereas (— logio P) 
around 0.43 if no signal is present, the value increases if a top is produced. In the j-j-j event 
class a (— logiQ P) of 2 is obtained for a cross-section atop of ~ 0.5 pb. This value can be 
compared with the 95% confidence level exclusion limit on the top production cross section 
at atop < 0.48 pb already derived by the HI experiment using the hadronic top decay channel 
only [30]. A deviation with three event classes with a (— log^Q P) > 3 would be found in this 
example for atop ~ 1.5 pb. 

The second test concerns the production of leptoquarks (LQs) [31]. Si/2,l and Vo,l type 
leptoquarks have been considered, which would mainly manifest themselves in the e-j and 
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j-u channels. A A coupling of 0.05 has been assumed and the sensitivity of the present analysis 
was tested for various LQ masses. The (— log^g P) values obtained from searches in the Mau 
distributions are summarised in figure |^ (bottom), for both the Si/2,l and Vq l LQ appearing in 
the e-j and e-j-j as well as the j-u and j-j-u channels, respectively. This analysis is sensitive to 
both types of leptoquarks up to masses of 240 — 250 GeV. These values can be compared with 
95% confidence level limits of 265 GeV for 5*1/2, l LQs and 240 GeV for Vo,l LQs, determined 
by dedicated analyses [32]. As for the case of single top production, the general search is thus 
found to have a sensitivity to leptoquark production which is comparable with that of dedicated 
searches. 
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UDject 


Energy Scale 


n 
O 


Identification 






(mrad) 


efficiency 


Jet 


2% 


5-10 




Electron 


0.7-3% 


1-3 


2-7% 


Photon 


0.7-3% 


1-3 


2-7% 


Muon 


5% 


3 


5% 



Table 1: Systematic uncertainties attributed to the measurement of energies and polar angles 
and to the identification efficiencies of particles. 



Process 



Uncertainty 



ep 
ep 
ep 
ep 
ep 
ep 
ep 



jjX and ep - 
jvX and ep - 
jjvX and ep 
ppX and ep - 
WX and ep - 
e'jX and ep - 
e'jp 



■jeX 

JjeX 
> eeX 
^WjX 

■ 



15% 

10% 

15% 

3% 

15% 

10% 

5% 



Table 2: Theoretical uncertainties attributed to the simulation of different SM processes. 





Mall 




event class 
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Nobs 


NsM 




SNsM 


P 


P 


Nobs 


NsM 


± 


SNsM 


P 


j-i 


0.38 


1 


0.035 




0.017 


0.036 


0.12 


1 


0.013 


± 


0.006 


0.013 




0.94 


111 


139 




21 


0.12 


0.021 


12 


31.2 


± 


5.1 


0.0028 


M-i 


0.67 


3 


1.07 


± 


0.25 


0.098 


0.29 


3 


0.70 


± 


0.23 


0.040 




0.34 


83 


116 


± 


14 


0.028 


0.22 


20 


36.7 


± 


6.2 


0.023 


e-u 


0.94 


5 


10.6 


± 


4.4 


0.17 


0.77 





2.1 


± 


0.8 


0.17 


e-e 


0.32 


3 


0.56 


± 


0.17 


0.023 


0.019 


3 


0.18 


± 


0.08 


0.0013 


e-n 


0.21 


4 


0.93 


± 


0.12 


0.016 


0.56 





2.6 


± 


0.5 


0.080 




0.069 


2 


0.14 


± 


0.04 


0.010 


0.036 


2 


0.11 


± 


0.03 


0.0060 


j-l 


0.52 


3 


10.8 


± 


3.7 


0.052 


0.77 





2.5 


± 


1.0 


0.13 


e-7 


0.38 


9 


19.2 


± 


2.0 


0.014 


0.64 


8 


15.7 


± 


1.9 


0.040 


7-7 


0.47 


1 


0.16 


± 


0.09 


0.15 


0.31 


1 


0.11 


± 


0.09 


0.12 


j-j-j 


0.41 


12 


5.9 


± 


2.0 


0.050 


0.58 


14 


7.8 


± 


2.5 


0.077 


e-j-j 


0.69 


39 


59.6 


± 


10.7 


0.058 


0.085 


9 


23.9 


± 


4.4 


0.0072 


3-3-^ 


0.62 


5 


1.79 




0.41 


0.043 


0.51 


5 


1.74 


± 


0.45 


0.040 


e-j-u 


0.090 


2 


0.19 




0.05 


0.016 


0.16 


2 


0.28 


± 


0.06 


0.034 


fl-j-U 


9.7- IQ-^ 


3 


0.19 


± 


0.05 


0.0011 


1.0 ■ 10"-^ 


3 


0.068 


± 


0.029 


7.5 ■ 10"^ 


3-3-1 


0.27 


1 


0.074 




0.048 


0.076 


0.36 


1 


0.15 


± 


0.10 


0.15 


e-i-7 


0.47 


1 


5.7 




1.6 


0.050 


0.39 


1 


5.6 


± 


1.4 


0.045 


e-j-j-j 


0.98 





1.6 




0.5 


0.23 


0.87 


1 


0.18 


± 


0.06 


0.17 


3-3-3-1^ 


0.33 


1 


0.084 




0.045 


0.083 


0.20 


2 


0.31 


± 


0.14 


0.044 



Table 3: The P values, the number of data events Nobs and the SM expectation Nsm for the 
region derived by the search algorithm using the Mau and ^ Pt distributions for event classes 
containing at least one event and taken into account in the statistical procedure. The p value in 
the selected region is also presented. 
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n 


Pi 


Pl.5 


P2 


P2.5 


P3 


P3.5 


Pa 


-P4.5 


P5 


1 


95% 


65% 


28% 


9% 


3% 


0.9% 


0.2% 


0.1% 


<0.05% 


2 


79% 


28% 


4% 


0.6% 


0.1% 


<0.05% 








3 


53% 


8% 


0.4% 


0.05% 


<0.05% 











Table 4: The probability PjJ to find at least n event classes with a — log P value greater than X. 
The values are applicable to both the M^n and ^ analyses. 
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Figure 1 : The data and the SM expectation for all event classes with a SM expectation greater 
than 0.01 events. The analysed data sample corresponds to an integrated luminosity of 1 17 pb^^. 
The error bands on the predictions include model uncertainties and experimental systematic 
errors added in quadrature. 
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H1 General Search - E Pj Distributions 
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Figure 2: The number of data events and the SM expectation as a function of ^ Pt for classes 
with at least one event. The shaded areas show the regions of greatest deviation chosen by the 
search algorithm. No search is performed for the j-j-j-j and e-j-j-j-j classes. 
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H1 General Search 
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Figure 3: The number of data events and the SM expectation as a function of Af^n for event 
classes with at least one event. The shaded areas show the regions of greatest deviation chosen 
by the search algorithm. No search is performed for the j-j-j-j and e-j-j-j-j classes. 
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Figure 4: The calorimetric energy deposits of the e-j-j-j-j event as a function of 7] and (f). 
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Figure 5: The — log^o P values for the data event classes and the expected distribution from 
MC experiments, as derived by investigating the Af^n distributions (top) and Pt distributions 
(bottom) with the search algorithm. 
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Sensitivity to Anomalous Top Production 
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Sensitivity to Leptoquark Production 
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Figure 6: The mean value of — log^Q P as derived from MC experiments which include a top 
signal with a cross section atop (top) and a LQ signal with a mass M^g and a A coupling equal 
to 0.05 (bottom), using the distributions of ^ Pt and Mau, respectively. The small arrows on 
the bottom figure indicate that these values should be treated as lower limits. 
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